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Green's  Function  Calculation  of  a Metal-Metal  Interface* 


A,  Yaniv 


Department  of  Physics,  University  of  California  San  Diego 
La  Jolla,  California  92093 


Abstract 

The  properties  of  metal-metal  interfaces  are  of  interest  in  many  metallurgical 
applications.  These  include  grain  boundaries,  crack  growth,  friction  and  bime- 
tallic adhesion.  The  present  work  is  a study  of  the  electronic  properties  of 
a simple  bimetallic  junction.  The  method  employed  to  investigate  the  inter- 
face are  the  Green's  function  technique  and  the  phase  shifts  method.  We  cal- 
culate the  Green's  function  of  a simple  junction  within  the  tight -binding 
approximation.  The  conditions  for  the  occurance  of  bound  states  is  deduced 
from  the  poles  of  the  Green's  function.  It  is  shown  that  there  are  three  types 
of  wave  functions  associated  with  the  interface.  The  first  extends  throughout 
the  entire  crystal,  the  second  extends  on  one  side  of  the  interface  only, 
whereas  the  third  kind  is  localized  near  the  interface.  Using  the  phase 
shifts  method  we  derive  an  expression  for  the  change  in  density  of  states  due 
to  the  creation  of  the  interface.  From  this  expression  we  derive  the  correspond- 
ing single-particle  contribution  to  the  interface  energy  and  the  interface  specific 
heat . 

* Supported  in  part  by  the  Office  of  Naval  Research  and  the  National  Science 
Foundation . 


I.  Introduction 


The  properties  of  metal-metal  interfaces,  are  of  interest  in  many  metalurgi- 
cal  applications.  These  include,  for  example,  grain  boundaries,  crack  growth, 
friction  and  bimetallic  adhesion.  The  general  metal-metal  interface  is  a much 
more  complicated  problem  than  the  corresponding  metal-vacuum  interface.  Recently 
there  has  been,  however,  some  success  in  applying  certain  methods  developed 
for  surface  phenomena,  to  the  investigation  of  bimetallic  interfaces. 

The  electronic  properties  of  1 -dimens iona 1 bimetallic  junctions  were  dis- 
cussed by  several  authors  using  various  types  of  models.  Aerts^^^  studied 
the  electronic  structure  of  a 1-dimensional  Kronig-Penney  model,  in  the  limit 
of  a 6-like  potential.  In  this  work  he  established  the  possibility  of  the 

existence  of  bound  interface  states.  Very  recently,  a relativistic  treatment 

(2 ) 

of  a similar  model  was  given  by  Kandilarov  and  Detcheva  . A tight-binding 

approach,  which  was  used  successfully  to  several  surface  problems,  was  also 

(3) 

applied  to  the  1 -d imens iona 1 interface.  Davison  and  Cheng  investigated  the 

electronic  properties  of  such  a system  using  the  molecular-orbita’  method. 

The  model  they  use  associates  a single  s-type  orbital  with  each  atom.  A 

(4) 

similar  model  was  studied  by  Allan  and  Lannoo  , using  an  approximate  Gaussian, 
density  of  states,  having  the  correct  second  moment.  (Green's  function  formalism 
was  used  by  Foo  and  Wong^^^  to  study  the  interface  states  of  a 1-dimensional 
sp -hybrid  junction. 

As  far  as  we  know,  the  study  of  three-dimensional  bimetallic  interfaces 
was  carried  out  by  the  density-functional  formalism  only.  This  formalism, 
developed  by  Hobenbery,  Kohn,  and  Sham^^^’^^\  was  applied  recently  to  the 

bimetallic  interfaces  by  Bennett  and  Duke^^^’^^\  Ferrante  and  Smith 

(12)  (13) 

Rouhani  and  Schuttler  , and  by  Mehrotra,  Pant  and  Das 

Whereas  the  density-functional  formalism  is  applicable  mainly  to  simple 
metals,  the  tight -binding  approximation  is  more  suitable  for  the  description 
of  transition  metals.  Since  the  purpose  of  the  present  work  is  to  investigate 
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tho  electronic  structure  of  a bimetallic  junction,  formed  from  two  transition 

metals,  wo  shall  use  the  tight-binding  approach.  The  model  we  consider  is  a 

highly  simplified  one.  The  two  metals,  on  each  side  of  the  junction,  are 

described  by  an  s-type  tight-binding  Hamiltonians.  The  interface  we  consider 

is  formed  by  bringing  together  two  semi-infinite  crystals,  and  creating  bonds 

between  the  atoms  on  the  two  surfaces.  We  assume  that  the  two  semi-infinite 

cyrstals  have  the  same,  two-dimensional,  translation  symmetry  parallel  to  the 

interface.  The  electronic  properties  of  this  model  are  investigated  by  using 

(14) 

the  Green's  function  method,  described  in  detail  by  Kalkstein  and  Soven 

The  details  of  the  tight-binding  model  are  outlined  in  Section  (2). 

Within  this  model  we  allow  for  a change  in  the  self-consistent  potential  of 
the  electrons  near  the  interface.  The  diagonal  matrix  elements  of  the  Green's 
function  are  calculated  in  Section  (3)  by  considering  the  formation  of  the 
interface  as  a perturbation  on  the  two  semi-infinite  crystals.  This  is  accom- 
plished by  the  use  of  Dyson's  equation.  Using  the  expression  for  the  Green's 
function  we  discuss  the  electronic  structure  of  the  interface  in  Section  (4). 
We  show  that  there  are  three  types  of  wave  functions  associated  with  the 
interface.  The  first  one  extends  throughout  the  entire  system,  the  second 
type  extends  on  one  side  of  the  junction  only,  and  the  third,  associated  with 
bound  states,  is  localized  near  the  interface.  The  behavior  of  the  bound 
states  as  a function  of  the  coupling  constant  between  the  two  metals  is  also 
discussed.  The  "phase  shifts"  method  is  applied  in  Section  (5)  to  determine 
the  change  in  the  total  density  of  states  due  to  the  creation  of  the  inter- 
face. Applying  this  result  we  obtain  an  expression  for  the  single-particle 
energy  contribution  to  the  interface  energy  and  the  corresponding  contribution 
to  the  electronic  specific  heat. 
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II  The  Model 


Consider  the  formation  of  a metal-metal  interface  by  bringing  together 
two,  semi-infinite,  metallic  crystals.  As  soon  as  a contact  is  formed,  elec- 
trons will  flow  from  the  metal  having  the  higher  Fermi  energy  to  the  one  having 
the  lower  energy.  This  flow  of  electrons  stops  when  the  potential  energy 
difference  between  the  two  sides  of  the  junction,  which  is  created  by  the  dipole 
layer  produced  at  the  interface,  is  equal  in  magnitude  and  opposite  in  sign 
to  the  difference  between  the  two  Fermi  levels.  Thus,  the  Fermi  levels  of  the 
two  metals  are  aligned,  due  to  the  interface  dipole-layer,  when  the  bimetallic 
junction  is  formed.  If  -Av  and  +Av  are  the  electrostatic  potential  created 
by  the  dipole  layer  on  the  right  and  the  left  hand  sides  far  away  from  the 
interface,  then  we  have: 

-eAv  = ^ (2.1) 

d b 

where  -e  is  the  charge  of  the  electron  and  E^  and  E^  are  the  original 
Fermi  levels  of  metals  a and  b , located  to  the  right  and  to  the  left  of 
the  junction,  respectively.  The  common  Fermi  level,  after  the  junction  is 
formed,  is  given  by: 

= E^^  + eAv 

= - eAv 

= % (E^^  + E^*")  (2.2) 

In  order  to  describe  the  electronic  properties  of  the  two  bulk  metals, 
the  two  metal-vacuum  interfaces,  and  the  metal-metal  interface,  we  apply  the 
t ight -binding  approximation  for  the  various  Hamiltonians.  For  simplicity  we 
associate  one  s-type  orbital  with  each  lattice  site.  It  is  also  assumed  that 
each  orbital  has  a g-fold  degeneracy  in  order  to  partially  account  for  the  10- 
fold  degeneracy  of  the  d orbitals  in  transition  metals.  Therefore,  the  Hamiltonians 
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I 


L 


of  the  two  bulk  metals,  and  , are  given  respectively  by 


ij 


a t 
t . . a .a  . 
ij  i J 


H,  - E t . . b . li 


i.J 


t 
.1 
I j 


(2.3) 


where  a . , a . and 
1 J 


b.^  b, 
1 ’ 


are  the  creation  and  destruction  operators  of  the 

orbitals,  associated  with  metals  a and  b localized  near  sites  i and  j. 

The  prime  on  the  summation  sign  denotes  a summation  over  nearest  neighbors 

3 h 

only.  The  matrix  elements  t..  and  t..  are  given  by 
^ ij  ij 

fE  if  i = j 
oa 


tT.  = • ilH  = 
ij  a ' 


or  E*^  if  i and  j are  nearest  neighbors 
0 otherwise 


(2.4a) 


E . if  i = j 
ob 

^Ib  ^Ib 


if  i and  j are  nearest  neighbors  (2.4b) 

lO  otherwise 

To  cleave  the  crystal  along  a given  low-order  crys ta lographic  plane,  we  have 

to  break  the  bonds  between  two  adjacent  planes,  parallel  to  the  corresponding 

direction.  We  neglect  any  geometrical  reconstruction  of  the  crystal  due  to  the 

3 b 

formation  of  the  surface  and  assume  that  the  transfer  integrals  t..  and  t.. 

ij  ij 

with  i j , have  the  same  value  as  in  the  bulk  systems,  provided  that  both 
sites  i and  j are  occupied.  In  addition,  we  assume  that  the  effect  of  the 
surface  on  the  redistribution  of  the  electronic  charge  near  the  surface  can  be 
described  as  a change  in  the  self-consistent  potential  of  the  electrons  near 
the  first  surface  layer  only.  Therefore,  the  diagonal  matrix  elements  of  sites 
located  on  the  surface  plane  will  be  different  from  the  corresponding  bulk  value. 
Thus,  for  the  metal-vacuum  systems  we  have 


i> 

= E 

+ U 

(2. 

5a) 

03 

3 

o 

1,0 

i> 

= E u 
ob 

1,-1 

(2. 

5b) 

where  0^  and  denote  the  change  in  the  self-consistent  potential  of  the 
electrons  near  the  surface.  In  the  above  expressions  we  assumed  that  i = 0 and 
i = -1  denote  surface  sites  of  metal  a and  b , respectively. 
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In  order  to  form  the  metal-metal  interface  we  start  from  two  semi-infinite  metal 
vacuum  interfaces,  having  the  same  translational  symmetry  parallel  to  the  sur- 
face, and  introduce  a coupling  between  the  two  surface  layers 


(2.6) 


where  R_^  and  ^ denote  two  nearest-neighbors  sites  located  on  the  two  sides 
of  the  interface.  We  also  note  that  the  dipole  layer  created  at  the  interface 
will  change  the  self-consistent  potentials  on  the  two  sides  of  the  junction. 

If  we  assume  that  charge  redistribution  is  confined  to  the  near  vicinity  of  the 
interface  we  can  write 


■-R . I H Ir.>  = E - eAV  + U 
a — 1 oa  a 


1,0 


(2.7a) 


(2.7b) 


where  i - 0 and  j - -1  denote  interface  sites  of  metals  a and  b , respec- 
tively. In  the  above  expressions  U and  U,  denote  the  change  in  the  self-con- 
sistent  potential  of  electrons  near  the  interface  relative  to  the  bulk. 
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III.  The  Bimetallic  Interface  Green's  Function 


As  we  already  noted  before,  we  assume  that  the  two  crystals  forming  the  inter- 
face have  the  same  translational  symmetry  parallel  to  the  interface.  As  a 
result  the  wave  vector  parallel  to  the  interface  is  a good  quantum  number 

in  such  a system.  In  the  following  we  assume  that  the  interface  is  in  the 
zy  plane,  and  that  metals  a and  b are  to  right  and  to  the  left  respectively. 
The  integers  m,  n,  ...  will  be  used  to  denote  the  various  planes  parallel  to 

the  interface.  Metal  a is  assumed  to  occupy  the  planes  m '•  0 whereas  m -1 

(14) 

planes  are  occupied  by  b atoms.  Following  Kalkstein  and  Soven  , we  denote 

by  _T^  the  translation  vector  parallel  to  the  interface  which  brings  the  atoms 
in  the  n-th  plane  to  coincide  with  those  on  the  plane  n = 0.  The  general 
coordinate  of  an  atom  on  the  n-th  plane  is  given  by: 


R = R 1 + R + T (3.1) 

n n II  n 


where  R j.  is  the  distance  between  the  planes  o and  n,  and  R is  a general 
n " 

translation  vector  parallel  to  the  interface.  Using  the  localized  orbitals 

|R  x’  + T > we  define  the  mixed  Bloch-Wannier  representation  by  the  follow- 
n ti 

ing  two-dimensional  Bloch  sum 


|k„,n>  = (N„)‘ 


R-. 


IV 


R-. 


+ T > e 


(3.2) 


where  N„  is  the  number  of  atoms  in  the  planes  parallel  to  the  interface.  These 
functions  are  localized  near  the  n-th  plane.  Since  the  translational  symmetry 
parallel  to  the  interface  is  preserved,  the  eigenfunctions  of  the  interface 
Hamiltonian  can  be  expanded  in  terms  of  the  Bloch-Wannier  functions  (3.2)  as 
follows 


k„>  = I 
n 


"nk„ 


Iki..  n> 


(3.3) 


where  j stands  for  any  quantum  number  needed  to  specify  the  eigenstates 
in  addition  to  k„. 


•i 


a 


{ 
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Using  the  expression  (3.3)  we  can  define  various  densities  of  states  for  the 
system  with  the  interface.  The  local  density  of  states,  near  the  n-th  plane, 
with  a definite  value  of  k„  is  given  in  terms  of  the  coefficients  c^j^  by 


(3.4) 


The  local  density  of  states  (per  atom)  near  the  n-th  plane  is 


0„(E) 


r ^ 

Nil  , 


(E,k„) 


= -L 

Nil 


j .Is... 


^nk„ 


(3.5) 


The  sum  over  ki,  in  this  expression  goes  over  all  values  of  k,i  in  the  two- 
dimensional  Brillouin  zone  defined  by  the  crystal  geometry  parallel  to  the 
interface.  Summing  the  expression  (3.5)  over  the  various  crystal  planes  we 
obtain  the  total  density  of  states.  The  density  of  states  per  atom  is  thus 
given  by 


0(E)  . 1-  £ (E)  . i E I cl  |S(E  - E ) (3.6) 

^ j »k|i 

n 

where  N is  the  number  of  planes  parallel  to  the  interface  and  N is  the  total 
number  of  atoms  given  by  N = N„  N 


In  order  to  derive  thevarious  densities  of  states  of  the  bimetallic  interface  we 
shall  use  the  Green's  function  technique.  Let  H and  be  the  Hamiltonians 

of  the  interfaced  crystal  and  of  the  free  metal-vacuum  interfaces  respectively. 
The  metal-metal  interface  Green's  operator  is  defined  by 


(E  - i6  - H)  G = 1 


(3.7) 


'where  E is  the  energy  and  6 is  a positive  infinitesimal.  The  surface  Green's 
operator  G°  is  related  to  the  Hamiltonian  H^  by  a similar  equation.  The 
Green's  operator  of  the  interface  is  related  to  the  surface  Green's  operator  via 
Dyson's  equation: 
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(3.8) 


G = G°  + G°VG 

where  V is  the  perturbation  necessary  to  create  the  metal-metal  interface 
from  the  free  surfaces  of  metals  a and  b , i.  e. 

V = H - H (3.9) 

o 

We  note  that  because  of  the  translational  symmetry  parallel  to  the  interface, 

G,  G°  and  V will  all  be  diagonal  in  the  wave  vector  index  k„  , in  the  Bloch- 
Wannier  representation  (3.2).  Omitting  the  corresponding  6-function, 

6 (k,,  - k„'),  we  use  the  notations  G (m,n;  k,,)  , G°  (m,n;  k„)  and  V (m,n;  k„) 
for  the  matrix  elements  of  G,  G°  and  V in  the  Bloch-Wannier  representation. 

To  further  simplify  the  notation,  we  shall  generally  omit  the  explicit  k,,  depen- 
dence in  these  expressions. 

The  various  densities  of  states  defined  earlier  by  Eqs  (3.4)  - (3.6),  are 
simply  related  to  the  imaginery  part  of  the  diagonal  matrix  elements  G (m,m;  k,,) 
by 

p(E)  ~ ImTr  G = ^ Im  Z G (m,m;  k„) 

m,k„ 

Pn(^)  “ G(n,n;  k„) 

li" 

Pj^(E,  k>t)  IraG  (n,n;  k„) 

We  turn  now  to  the  evaluation  of  the  interface  Green's  function  from  Dyson's 
equation  (3.8).  This  operator  equation  is  reduced  to  the  following  algebraic 
equation  in  the  Bloch-Wannier  representation: 

G(m,n)  =G°(m,n)  +S  G°  (m,  £)V  (.^’,r)G  (r  ,n)  (3.13) 

i,r 

From  the  discussion  of  the  model  given  in  Section  (2),  it  is  obvious  that  the 
perturbation  potential  V has  off-diagonal  matrix  elements,  which  couple  the 


(3.10) 

(3.11) 

(3.12) 


-8- 


3 


i 


i' 


interface  planes  n = 0 and  n = -1.  In  addition,  on  each  side  of  the  inter- 
face V has  two  types  of  diagonal  matrix  elements.  The  first  one,  due  to  the 
interface  dipole  layer,  is  given  by  iteAv.  This  is  just  the  perturbation 
necessary  to  align  the  Fermi  levels  on  the  two  sides  of  the  junction,  Eq . (2.1). 
The  second  type  of  diagonal  matrix  elements  is  due  to  the  change  in  the  self- 
consistent  potential  of  the  electrons  near  the  interface.  This  is  given  by 

U - U and  U - U,  where  U , U and  U U,  are  the  self-consistent  potential 
aabb  a’ab’b 

relative  to  the  bulk  of  interface  and  surface  electrons  of  metals  a and  b, 
respectively . 

As  we  shall  see  later  on,  it  is  possible  to  determine  the  surface  Green's  func- 
tions for  an  arbitrary  value  of  the  surface  self-consistent  potentials  and 

. It  is , therefore,  simpler  to  start  from  a fictitious,  intermediate  surface 
problem  where  the  surface  self-consistent  potentials  and  U^,  have  already 

the  correct  interface  values,  u and  U,  , respectively.  Using  this  sytem 
as  our  starting  point,  the  perturbation  potential  necessary  to  create  the  metal- 
metal  interface  is  simpler  than  in  the  original  problem.  It  consists  of  diagonal 
terms  teAv,  which  align  the  Fermi  levels  of  the  two  metals  and  of  off-diagonal 
terms,  which  couple  the  two  semi-infinite  crystals.  Thus,  the  only  non-vanishing 
matrix  elements  of  the  perturbation  are  given  by 

!-eAv  for  m - 0 

(3.14) 

-heAv  for  m s -1 


and 


V(-1,0)  = V(0,-D*  = g(k„)  (3.15) 

where  3 and  cp  are  derived  from  the  relation 

V(-1,0)  = Z'  <R  , |v|r  > e^-"  ■ 

R 

-o 

In  this  expression,  is  a translation  vector  parallel  to  the  interface  which 

brings  the  atomic  sites  of  the  two  interface  planes  to  coincide. 
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We  start  now  from  the  intermediate  surface  problem  discussed  above, whose  Green's 
function  we  denote  by  G°,  and  apply  the  perturbation  in  two  steps.  First  we 
align  the  Fermi  levels  of  the  two  metals  by  applying  a constant  electrostatic 
potentials  -Av  on  metal  a and  +AV  on  metal  b.  The  only  effect  of  applying 
these  constant  potentials  on  the  surface  green's  functions  is  a shift  in  the 
corresponding  energies.  Explicitly  we  have 


G°  (E)  = G°  (E-heAv) 

3 3 


G°  (E)  = G°  (E-eAv) 


(3.17a) 


(3.17b) 


where  G denotesthe  value  of  the  Green's  functions  in  the  intermediate  problem, 
after  the  application  of  the  potentials  xAv.  To  simplify  the  notation  we  drop 
the  wiggle  from  these  Green's  functions  remembering  that  the  energies  have  to  be 
shifted  according  to  (3.17). 

At  this  stage  we  turn  on  the  second  part  of  the  perturbation  i.e.  the  off- 
diagonal  coupling  between  the  two  metals.  Due  to  the  localized  nature  of  this 
perturbation,  the  Dyson's  equation  (3.13)  is  greatly  simplified  and  we  obtain 


G(m,n)  =G°(m,n)  +G°(m,-1)  V(-1,0^  G(0,n)  +G°(m,0)  V (0, -1)G (-1  ,n)  (3.18) 

Using  the  fact  that  the  surface  Green's  function  G°(m,n)  vanishes  if  m and 
n refer  to  planes  on  opposite  sides  of  the  interface,  we  can  easily  solve  equa- 
tion (3.18)  for  the  perturbed  Green's  function.  In  this  way  we  get 

G(m,n)  =G°(m,n)  +G°(m,0)  [ V (0 , -1)G°  (-1  ,n)  + | V (0 , -1)  | V (-1 , -1)  x 

G°(0,n)]  [l  - 1V(0,-1)|^  G°(0,0)G°(-l,-l)j 


for  m '0 


(3.19a) 
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(3.19b) 


for  m s -1, 


G(m,n)  = C°(m,n)  + (m, -1)  [ V (-1 ,0)  G°(0,n)  + | V (0, -1)  | ^G° (0,0)  x 

G°(-l,n)|  [l-  |V(0,-1)!^  G°(0,0)  G°(-l,-l)  ] 


In  particular,  the  diagonal  matrix  elements  of  the  interface  Green's  function 
are  given  by 

G(m,m)  =G°(m,m)  +B^G°(-1,-1)  G° (m,0)G° (0,m)  [ 1-B V (0 ,0)G°  (-1, -1 )] 
for  m = 0 (3 


G(m,m)  = G°(m,m)  + B^G°(o,o)  G°(m,-1)  G°(-l,m)  | 1-B^  G°  (0, 0)G°  (-1 , -1)  | 

for  m s -1  (3 .201 

In  the  above  expressions  B is  related  to  the  off-diagonal  coupling  between 
the  two  metals  through  Eq.  (3.15),  i.e. 

3^  = |v(0,-l)|^ 


The  diagonal  matrix  elements  of  the  surface  Green's  function  were  calculated 
by  Kalkstein  and  Soven^^^^  Using  the  same  method  one  can  also  derive  the 
general  matrix  elements  of  this  Green's  function.  In  Appendix  A we  show  that 
these  matrix  elements  are  given  by 


G (m,n)  = iti, 


■1  -i(m-n)9  1 / 

• ( 


a 3 


\ 2T  / ip,- 


iU+(uj  -2U  ) ) 

a a a ( , 


p-(iD  -2U  ) 
'^a  a a 


for  m,  ni  0 


tm-n  1 


ml  + nl -2 


G (m,n)  = iu-L  e 


for  m,n  s -1 


and 

G°(m,n)  = 0 otherwise.  (3.21c) 

In  these  expressions  p,  and  o)  are  defined  by 
(4T^-iu^)^  for  s 4t^ 
isign(uJ)  ((D^-4T^)^  for  ^ > 4T^ 

and 

ID  = E - W(k„) 

where  W,  T,  and  6 are  related  to  the  matrix  elements  of  the  bulk  Hamiltonians 
of  the  two  metals  in  the  Bloch-Wannier  representation  as  follows 

”a(b)  = ■ nk„I  (bulk)  |nk„> 

Tg(b)(!iM)  e^®a(b)^-"^  = >^nk„|H^^j^^  (bulk)  \ n+l,k„> 

In  order  to  apply  the  expressions  (3.21)  to  the  solution  of  the  interface 
problem,  we  replace  first  and  by  the  corresponding  interface  values 

and  Uj^,  respectively.  In  addition  we  have  to  shift  the  energies  according 
to  Eq  . (3.17)  in  order  to  align  the  Fermi  levels  of  the  two  metals.  This  can 
be  achieved  by  redefining  the  id-s  in  the  following  way 

CD  = E - W (k„)  + eAv  (3.22a) 

3 3 

cDj^  = E - Wj^(k„)  - eAv  (3.22b) 

Substituting  the  explicit  expression  (3.21)  for  the  surface  Green's  function 
into  Eq.  (3.20),  we  obtain  that  the  diagonal  matrix  elements  of  the  interface 
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Green's  function  are  given  by 


We  recall  that  and  in  these  equations  are  defined  by  the  energy-shifted 

expressions  (3.22a)  - (3.22b),  respectively. 
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[V.  The  Electronic  Structure  of  the  Interface 


In  this  section  we  apply  the  Green's  function  derived  in  the  preceding  section 
in  order  to  investigate  the  electronic  structure  of  the  interface.  From  the 
explicit  expression  (3.23)  it  follows  that  the  diagonal  matrix  elements  of 
the  interface  Green's  function  have  a non-vanishing  , continuous  imaginery  part 
for  energies  which  are  either  inside  the  band  of  metal  a or  inside  the  band 
of  metal  b (shifted,  of  course,  according  to  (3.17)).  Therefore,  as  follows 
from  Eq . (3.12)  the  band  width  of  the  crystal  with  the  interface  is  the  union 
of  the  shifted  bands  of  the  two  separate  crystals.  However,  the  wave  functions 
of  the  combined  crystal  can  be  classified  into  three  distinct  classes  according 
to  their  different  localization  properties.  To  facilitate  our  further  discussion 
we  define  the  k„  sub-bands  of  the  two  metals  as  the  band  structure  obtained  by 
the  intersection  of  E (k)  and  E (k)  with  the  plane  k,,  = constant,  where  E (k) 

3D  3 

and  (ji)  the  single-particle  energy  spectra  of  the  two  bulk  crystals. 

The  first  type  of  states  have  wave  functions  which  extend  throughout  the  entire 
crystal.  As  can  be  seen  from  Eqs . (3.4)  and  (3.23),  this  kind  of  behavior  is 
associated  with  states  whose  energy  lies  in  the  k„  sub-band  of  the  two  metals, 
i.e.  their  energy  and  wave  vector  satisfy  the  relations 

Ie  - W (k„)  + eAv|  <2T  (k„) 

a ““  a 

1e  - Wj^(k„)  - eavl  <2Tj^(k„) 

The  second  type  of  states  have  wave  functions  which  extend  to  ^ on  only  one 
side  of  the  interface,  and  which  decay  exponentially  with  the  distance  from 
the  interface  on  the  other  side.  Using  the  expression  (3.23)  for  the  Green's 
function  it  is  not  hard  to  see  that  this  behavior  occurs  when  the  energy  lies  in 
the  k„  sub-band  of  one  metal  but  outside  the  corresponding  sub-band  of  the  other. 
Thus,  for  energies  and  wave  vectors  that  satisfy 

(4.2) 


j 

(4.1) 


|e  - W (k„)  + eAv|  2T  (k„) 

u a 

|e  - Wj^(k„)  - eAvl  >2T^(k„) 
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i 
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the  corresponding  wave  functions  are  Bloch-like  inside  metal  a but  decay 
exponentially  from  the  interface  inside  metal  b.  In  a similar  way,  the  wave 
functions  of  states  whose  energies  and  wave  vectors  satisfy 

|e  - W (k„)  + eAv|>2T  (k„) 

cl  3 ~ 

I I 

|e  - W|^(k„)  - eAv|<2T^(k„) 

extend  to  infinity  on  the  b side  of  the  interface,  but  decay  exponentially 
on  the  a side.  The  decay  coefficient  of  the  wave  function  is  determined  by 
die  energy  measured  relative  to  the  corresponding  k,,  sub-band  center.  If 
we  express  this  energy  in  units  of  the  sub-band  half-width,  and  write 

ei)  (k„)  = 02T(k„)  , \a\'A  (4.4) 


where  ^(k„)  and  T (k„)  refer  to  the  corresponding  values  on  the  side  of  the 
interface  where  the  wave  function  decays  exponentially,  we  obtain  from  Eqs . 
(3.4)  and  (3.23)  that  on  the  decaying  side  the  wave  function  on  plane  m is 
proportionaly  to  e . The  decay  coefficient  \>0  is  given  by 


X = - ^n  O'  - s ign  (O') 


Thus,  the  further  the  energy  is  from  the  center  of  the  sub-band,  the  stronger 
is  the  exponential  decay  of  the  corresponding  wave  function. 


We  note  that  the  behavior  of  the  wave  functions  extending  throughout  the  entire 
system  and  those  extending  only  on  one  side  of  the  interface  is  as  expected. 

If  we  try  to  propagate  a wave  through  the  interfaced  crystal  this  wave  can 
propagate  from  one  side  to  the  other  only  if  its  frequency  is  in  the  common 
sub-bands.  If,  however,  the  frequency  is  in  the  sub-band  of  one  of  the  crystals 
but  outside  the  sub-band  of  the  other,  this  wave  cannot  penetrate  into  the 
second  crystal,  and  its  amplitude  will  decay  exponentially. 
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The  third  type  of  wave  functions  is  associated  with  the  existence  of  bound  inter- 
face states.  For  a given  k„,  the  energy  of  the  possible  bound  states  are 
determined  by  the  poles  of  the  Green's  function,  which  lie  outside  the  k,, 
sub-bands  of  the  two  metals  forming  the  junction.  Using  the  explicit  expression 
(3.23)  for  the  interface  Green's  function,  we  see  that  the  bound  states  ener- 
gies are  given  by  the  roots  of  the  equation: 

43^  - |^sign(uj^)|Ig+  (iJg  - 2uJign(cUj^)^^+  - 2uJ  = 0 (4.6) 

outside  the  k,,  sub-bands.  In  this  equation  |i.  is  defined  by 
= (u)2  . 4T^)^ 


As  can  be  seen  from  Eq.  (3.23),  the  wave  function  of  an  interface  bound  state 
is  localized  near  the  interface  and  decays  exponentially  with  distance  on  both 
sides  of  the  interface , Thus,  the  bound  electron  is  free  to  propagate  parallel 
to  the  interface,  but  is  confined  to  a finite  stripe  of  width  >,^  +>^  (where 
>.  and  X,  are  given  by  Eq.  (4.9))  in  the  direction  perpendicular  to  the  inter- 
face.  The  decay  coefficients  of  the  bound  state  wave  function  will  be  different, 
in  general,  on  the  two  sides  of  the  interface.  If  we  express  the  bound  state 
energy,  relative  to  the  center  of  the  two  k,,  sub-bands,  in  terms  of  the  corres- 
ponding sub-bands  half-widths,  we  can  write 


0)  = a 2T 


a 2T 
b b’ 


(4.7) 


where 


- E°  - W + eAv  ; - E°  - W,  - eAv  (4.8) 

a a b b 

and  E°  is  the  energy  of  the  bound  state.  It  is  easy  to  see  from  Eq . (3.23) 

that  the  exponential  decay  coefficients  on  side  a and  side  b , \ and 

a b 

respectively,  are  given  by 


16- 


- sign(aj^)v4^2_j  j 


(4.9a) 


(4.9b) 


Thus,  the  greater  the  distance  of  the  bound  state  energy  from  the  center  of  the 
k„  sub-bands,  the  more  localized  is  the  corresponding  wave  function. 
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The  analytical  solution  of  Eq . (4.6)  for  the  bound  states  energies  is  not 
possible  in  general.  However,  by  a suitable  graphical  analysis,  one  can  deter- 
mine the  conditions  for  the  occurance  of  bound  states,  their  number  and  their 
position  with  respect  to  the  bands  as  a function  of  the  strength  of  the  coupling 
constant.  It  turns  out  to  be  very  convenient  to  analyze  the  interface  bound 
states  in  terms  of  the  properties  of  the  bound  surface  states  of  the  intermediate 
surface  problem  considered  in  the  preceeding  section  (i.e.  the  surface  problem 
having  a surface  self-consistent  potential  equal  to  the  respective  value  of 

the  interface  problem).  It  is  well  known  that  the  bound  surface  states  are 

(14) 

given  by  the  roots  of  the  equation  : 
sign(u))p,  + UJ-2U  = 0 


The  corresponding  surface  bound  state  energies  are  given  by 

o ^a^ 

Ea  = - eAv  + f 

n 

E = W (k„)  + eAv  + U + 

b b - U. 


(4.10a) 

(4.10b) 


provided  U >T  and  U,  I > T,  . For  U T there  are  no  surface  states. 

' a a ' b b 

In  the  following  analysis  we  shall  use  the  indices  1 and  2 to  denote  the  metal 
having  the  higher  sub-band  top  edge  and  the  lower  sub-band  bottom  edge,  respec- 
tively. 

The  behavior  of  the  interface  bound  states  can  be  described  as  follows  ; 


In  the  intermediate  surface  problem  there  are  no  bound  states  (i.e. 


I U KT  and 
a ' a 


Uj^i'-'T^)  there  will  be  no  bound  Interface  states  for  small 


values  of  the  coupling  constant  8 . As  the  coupling  constant  is  increased, 
a bound  state  will  appear  as  soon  as  the  critical  value  3^ is  reached,  where 


bJ  - ..in  i <Ti-»i)£2«Ux)] 


(4.11) 
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•w- 


and 


f . (E)-sign (oj.  )u  +UJ.-2U.  (4.12) 

1 1 1 1 

^ 2 

E max  and  E min  are  the  top  edge  and  the  bottom  edge  energies  of  the  higher 

2 

and  the  lower  of  the  two  k,,  sub-bands  respectively.  If  - ) f . (E  min)  < 

1 4/1 

(T,  - U-)f«(L  ) the  bound  state  will  appear  below  the  sub-bands.  If,  how- 

ever,  - (T2+U2)f  (E  min)  > (Tj^-Uj^)f2  (E  max)  the  bound  state  will  appear  above 
the  sub-bands.  The  bound  state  described  above  exists  as  long  as  p . When 
the  coupling  constant  is  further  increased  and  the  value  is  reached,  a 
second  bound  state  appears  where 

for  5 ^ 32  there  are  always  two  bound  states,  one  below  and  the  other  above 
the  sub-bands. 

b.  Suppose  that  in  the  intermediate  surface  problem  one  of  the  following 
four  situations  holds: 

1.  There  is  one  bound  state  above  the  sub-bands  and  the  two  sub-bands 
overlap. 

2.  There  is  one  bound  state  above  the  sub-bands,  which  is  associated 
with  the  metal  having  the  higher  sub-band  and  there  is  a gap  between  the 
two  sub-bands. 

3.  There  are  two  surface  states,  one  above  the  sub-bands  and  the  other 
lies  in  the  sub-band  of  the  second  metal  and  tiere  is  an  overlap  between 
the  sub -bands. 

4.  There  are  two  bound  states,  a gap  exists  between  the  sub-bands  and  the 
bound  state  of  the  lower  sub-band  lies  above  the  sub -bands,  whereas , the  one 

associated  with  the  upper  sub-band  falls  in  the  lower  sub-band. 

Then  the  behavior  of  the  interface  bound  states  is  as  follows:  For  small  coup- 

ling constant  there  will  be  one  bound  state  above  the  sub-bands.  This  state 
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develops  in  a continuous  way  from  the  corresponding  surface  state.  If 
is  the  energy  of  the  surface  state,  which  lies  above  the  sub-bands,  then  for 
small  3 the  energy  of  the  corresponding  interface  state  is  given  approxi- 
mately by 


43" 


f.(Ej) 


(1 


(4.14) 


where  fj  (E)  is  the  function  defined  by  (4.12).  As  long  as  3 is  in  the 
range  o^3-  0^,  where  P^is  defined  by 


83  = (T^+h^lfjCE^min)  (4.15) 

there  is  one  bound  state  above  the  sub-bands.  For  3‘-P3  there  are  two  bound 
interface  states,  one  above  and  the  other  below  the  sub -bands. 

c.  Suppose  that  in  the  intermediate  surface  problem  one  of  the  following 
four  situations  exists: 

1.  There  is  an  overlap  between  the  sub-bands,  a bound  state  exists  below 
the  sub-bands. 

2.  There  is  a gap  between  the  sub-bands,  there  is  one  surface  state 
below  the  sub-bands  which  is  associated  with  the  lower  sub-band, 

3.  There  is  an  overlap  between  the  sub-bands  and  there  are  two  surface 
states,  one  below  the  sub-bands  and  the  second  lies  inside  the  sub-band 
of  the  other  metal. 

4.  There  is  a gap  between  the  sub-bands  and  there  are  two  surface  states, 
the  one  associated  with  the  higher  sub-band  lies  below  the  sub-bands,  and 
the  other  falls  inside  the  higher  sub-band. 

Then  the  behavior  of  the  interface  bound  states  is  as  follows  : For  small 

3 there  will  be  one  bound  state  below  the  sub-bands,  with  an  energy  given 

approximately  by  (4.14).  When  the  coupling  constant  is  increased,  a critical 

value  0 , is  reached,  where  a second  bound  state  appears  above  the  sub-bands, 

4 

0 is  defined  by 
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(4.  If)) 


“4  - 4 

For  3‘ 3^  there  are,  thus,  two  bound  states,  one  below  and  the  other  above  the 
sub -bands . 


d.  Suppose  that  in  the  Intermediate  surface  problem  one  of  the  following 
situations  holds: 

1.  There  is  a gap  between  the  sub-bands  and  there  exists  one  surface 
state  above  the  sub-bands,  which  is  associated  with  the  lower  sub-band. 

2.  There  is  a gap  between  the  sub-bands  and  there  exist  two  surface 
states.  The  one  associated  with  the  higher  sub-band  is  located  above  the 
sub-bands  and  the  other  lies  inside  the  higher  sub-band. 

Then  the  behavior  of  the  interface  states  is  as  follows:  for  a small  coupling 
constant  there  is  one  interface  state  above  the  sub-bands.  When  S reaches 
the  critical  value  where 

“5  ■ "1  '2  <'■• 

(E^^in  being  the  bottom  edge  energy  of  the  higher  sub-band)  a second  bound 

state  appears  inside  the  gap.  This  bound  state  exists  as  long  as  the  coupling 

constant  is  in  the  range:  0s  3 s 6, , where  6 is  given  by 

5 6 6 

i (h-  “2>  h «Lx> 

2 

E is  the  energy  at  the  top  of  the  lower  sub-band.  For  3 > F,  the  bound 

max  6 

state  disappears  from  the  gap.  When  B reaches  the  value  , given  by  Eq. 
(4.15),  a new  interface  state  appears  below  the  sub-bands.  For  B > there 

are  two  interface  states,  one  below  and  the  other  above  the  sub-bands. 


e.  Suppose  that  in  the  intermediate  surface  problem  one  of  the  following 
situations  holds: 
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Then  the  behavior  of  the  interface  states  is  as  follows:  For  small  values  of 

the  coupling  constant  there  is  one  interface  state  below  the  sub-bands  whose 

energy  is  given  approximately  by  (4.14).  When  S is  increased,  and  reaches 

the  value  3 given  by  (4.18),  a new  interface  state  appears  in  the  gap.  This 

6 

bound  state  exists  in  the  gap  for  3^3^  where  is  given  by  (4.17). 

For  B ■ 8^  the  bound  state  disappears  from  the  gap  and  a new  bound  interface 


state  appears  above  the  sub-bands  for 


■ 5^,  where  is  given  by  (4.16). 


For  5 ■ ? there  are  two  bound  states,  one  below  and  the  other  above  the  sub- 

4 

bands . 


f.  Suppose  that  one  of  the  following  situations  holds  in  the  intermediate 
^ surface  problem; 

I 1.  There  is  a gap  between  the  sub-bands  and  there  is  a single  surface 

j state  which  lies  in  the  gap. 

2.  There  is  a gap  between  the  sub-bands  and  there  are  two  surface  states, 

one  in  the  gap  and  the  other  lies  inside  the  sub-band  of  the  other  metal. 


Then  the  behavior  of  the  interface  bound  state  is  as  follows:  For  small  values 

of  the  coupling  constant  there  is  a single  interface  state  which  lies  in  the 
gap.  When  3 is  increased  and  reaches  the^^lue  where 


3 = max 

7 


(T^-nj^) 


(E 


min 


).  ^ 


(T2-U2) 


(E 


max) 


(4.19) 


this  bound  state  disappears.  For  3 3^ there  appears  a new  interface  state 

below  the  sub-bands  where  is  given  by  (4.15).  If  S ' 3^  there  is  also  a 

bound  state  above  the  sub-bands,  where  3^  is  determined  by  (4,16). 
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R.  If  in  the  intermedi.ite  surface  prohlem  there  are  two  surface  states  below 

the  sub-bands,  and  a gap  exists  between  those  sub-bands,  then  for  small  values 

of  there  will  be  two  interface  states  below  the  sub-bands  (with  energies 

given  approximately  by  (4.14)).  When  reaches  the  value  given  by 

(4.15),  the  interface  state  having  the  higher  energy  disappears.  A new 

bound  interface  state  appears  in  the  gap  when  p reaches  the  value  6.  where 

0 

is  given  by  (4.18).  This  bound  state  exists  as  long  as  where 

o 6 5 

is  determined  by  (4.17).  For  S ■ the  interface  state  disappears  from 
the  gap  and  a new  bound  state  appears  above  the  sub-bands  when  S 5^,  where 
is  given  by  (4.16). 

h.  If  the  situation  in  the  intermediate  surface  problem  is  the  same  as  des- 

cribed in  (g.)  except  that  there  is  an  overlap  between  the  sub-bands,  there 
will  be  two  interface  states  below  the  sub-bands  for  0 P ' 6^ , where  5^ 

is  given  by  (4.15).  A new  interface  state  appears  above  the  sub-bands  for 

P ■ , where  5^  is  given  by  (4.16). 

i.  If  in  the  intermediate  surface  problem  there  are  two  surface  states 

above  the  sub-bands,  and  a gap  exists,  then  for  small  values  of  P there 

are  two  interface  states  above  the  sub-bands.  When  9 reaches  the  value  9, 

4 

given  by  (4.16),  the  interface  state  having  the  lower  energy  disappears. 

A new  bound  state  appears  inside  the  gap  when  9 reaches  the  value  given 

by  (4.17).  This  bound  interface  state  exists  as  long  as  ■■  9 ' where 

is  given  by  (4.18).  If  P is  further  increased,  the  bound  state  dis- 
appears from  the  gap  and  a new  interface  state  appears  below  the  sub-bands 
when  P 9^  where  9^  is  given  by  (4.15). 

j.  If  the  situation  in  the  intermediate  surface  problem  is  the  same  as 

described  in  (i)  except  that  there  is  no  gap  between  the  sub-bands,  there 

will  be  two  interface  states  above  the  sub-bands  for  0 s-  ® P,  where  8, 

4 4 

is  given  by  (4.16).  For  8 > 8^  the  interface  state  having  the  lower  energy 
disappears  and  a new  one  appears  below  the  sub-bands  when  ° > 9^  where 
9^  is  given  by  (4.15). 
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k.  If  in  the  intermediate  surface  problem  there  are  two  surface  states 
inside  the  gap,  then  for  small  values  of  the  coupling  constant  there 
will  be  two  interface  states  inside  the  gap.  Their  energies  are  given 
approximately  by  (4.14).  When  P reaches  the  value  where 

O 

el  - min  [i  (Ij-up  (Tj  + Uj)  (4.20) 


one  of  these  bound  states  disappears.  When  B is  further  increased  and 
reaches  the  value  p^  given  by 


Pg  - max 


k(r^ 


■V 


(E  ), 

max 


(T, 


Ui) 


(E  . ) 
min 


(4.21) 


The  other  bound  states  also  disappears  from  the  gap.  For  P p^,  where 
p^  is  given  by  (4.15),  an  interface  state  appears  below  the  sub-bands, 
whereas  for  P • , where  P^  is  given  by  (4.16),  there  is  a bound  state 

above  the  sub-bands. 


1.  If  in  the  intermediate  surface  problem  there  are  two  surface  states, 
one  below  and  the  other  above  the  sub-bands,  then  for  any  value  of  the 
coupling  constant  p there  will  be  two  interface  states,  one  below  and 
the  other  above  the  sub-bands. 


m.  If  in  the  intermediate  surface  problem  there  are  two  surface  states, 
one  in  the  gap  and  the  other  below  (above)  the  sub-bands,  then  for  small 
values  of  P there  will  also  be  two  interface  states,  one  in  the  gap 
and  the  other  below  (above)  the  sub-bands.  The  bound  state  in  the  gap 
disappears  for  P > P^  , where  By  is  given  by  (4.10).  For  P ■ p^ 

(P  '■  P^)  a new  bound  state  appears  above  (below)  the  sub-bands,  where 
P^  (P3)  is  given  by  (4.16)  ((4.15)). 

This  analysis  covers  the  various  possible  bound  interface  states. 


We  note  that  a general  property  of  the  bound  interface  states,  derived  from  Eq. 
(4.6),  is  that  whenever  two  interface  states  exist  simultaneously,  an  increase 
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in  the  coupling  constant  will  increase  the  energy  of  the  interface  state  having 
the  higher  energy  and  will  decrease  the  energy  of  the  lower  energy  bound  state. 

Another  remark  that  should  be  added  here  is  that  for  very  strong  coupling  between 
the  two  metals  there  are  always  two  interface  states,  one  below  and  the  other 
above  the  sub-bands.  The  asymptotic  energies  of  these  bound  states  are  given 
by 


E,  = U + W - eAv  + 5 
1 a a 


E_  - U,  + W + eAv  - 0 
z b b 


if  I ^ I Ujj  ^b  ^ 


E,  “ U,  + W + eAv  + 0 
1 D b 


E-  = U + W -eAv  - 0 
2 a a 


if  lu.  + W,  + eAv  1 > 1 U + W - eAv  > 
b b ' a a 
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I If  many-body  effects  are  neglected  the  total  energy  of  the  crystal  with  the 

i interface  is  just  the  sum  of  the  occupied  single-particle  energy  levels.  This, 

in  turn,  can  be  expressed  as  a corresponding  integral  over  the  system's  density 
of  states.  The  total  density  of  states  of  the  interfaced  crystal  can  be  obtained 
from  the  Green's  function  derived  in  Section  (3),  by  summing  the  imaginery  part 
over  the  various  planes  parallel  to  the  interface.  However,  in  order  to  deter- 
mine the  energy  needed  to  break  the  metal-metal  interface  into  two  metal-vacuum 
interfaces  we  need  to  know  the  change  in  the  total  density  of  states  due  to 
this  cleavage  process.  This  can  be  done  directly,  without  calculating  the 
density  of  states  of  the  two  systems  by  using  the  phase  shifts  method  described 
by  De  Witt^^^\  Callaway  and  Toulouse  This  method  can  be  summarized 

as  follows;  Let  H=H  + V be  the  perturbed  Hamiltonian,  and  let  E.  and  e.  be 

o ’ J J 

j the  eigenvalues  of  H and  H^  respectively.  The  change  in  density  of  states 

due  to  the  perturbation  V can  be  written  as 


I 


Ap(E)  = p (E)  - (E)  = Z 

j 


6(E- 


6(E-^^)j 


where  and  p are  the  unperturbed  and  the  perturbed  densities  of  states. 

This  can  be  also  be  expressed  in  the  following  way: 


A Art  \ Xth  d 

^ 


-Cn 


E-i6-E. 
L 

E-i6-e . 

J 


where  G and  G°  are  the  Green's  operators  of  the  perturbed  and  the  unperturbed 
Hamiltonians,  respectively.  Using  Dyson's  equation  we  can  express  G~^  in  terms 
of  G°.  In  this  way  we  get 


Ao  (E) 


Im  Tr  [^/n  (1-G°V)  j 
^ (1-G°V)] 


(3.1) 
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In  the  last  step  we  used  the  identity  Tr  j/nAj  = In  |^det  aJ  , where  A is  a 
general  operator.  The  application  of  i^q . (5.1)  is  especially  useful  for  a 
localized  perturbation,  when  det  (1-(’.°V)  can  be  expressed  as  a finite-order 
determinant . 


Let  us  turn  now  to  the  specific  problem  of  the  interface.  In  this  case  we 
start  from  the  two  semi-infinite  surface  systems  (the  a and  the  b surfaces), 
and  apply  the  perturbation  necessary  to  create  the  interface.  This  perturbation 
was  described  in  detail  in  Section  (2).  We  introduce  the  perturbation  in  two 
steps.  In  the  first  one  we  apply  a constant  electrostatic  potential  -Av  on 
metal  a and  +Av  on  metal  b . As  we  have  seen  before  the  application  of 
these  potentials  causes  a shift  in  the  densities  of  states  of  the  two  surface 
systems,  which  aligns  the  Fermi  levels  of  the  two  metals.  Since  each  metal  is 
electrically  neutral,  there  will  be  no  net  change  in  their  energies  due  to  the 
application  of  this  perturbation.  Our  second  and  final  step  is  to  apply  the 
remaining  perturbation  V needed  to  form  the  interface.  As  we  have  seen  in 
Sections  (2)  and  (3),  the  only  non-vanishing  matrix  elements  of  this  perturba- 

tion are  given  by 

V (0,0)  = U - 0^  (5.2a) 

a a 

V (-1,-1)  = Ub  “ ^b 


and 


V (-1,0)  = V(0,-1)*  = 3e^^  (5,2c) 

where,  as  before,  U , U.  , and  U , U,  are  the  self-consistent  potentials  near 
’ ’ a’  b’  a’  b < 

the  interface  and  near  the  surface,  respectively.  Thus,  in  the  present  case 
we  have 


det  (1-G°V) 


1 - G°(-l,-l)  V (-1,-1) 
-G°(0,0)  V (0-1) 


-G°(-l,-l)  V (-1,0) 
1 - G°(0,0)  V (0,0) 


(5.3) 
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where  G°(0,0)  and  are  the  corresponding  surface  Green's  functions 

of  metals  a and  b respectively  with  the  energies  shifted  according  to 
(3,17).  We  note  that  the  surface  self-consistent  potentials  in  G°(0,0)  and  G°  (-1,-1) 
are  0^  and  respectively  and  not  and  as  was  the  case  in  the  preceding 
section.  Since  G°  and  V are  diagonal  in  the  wave  vector  k„  , the  deter- 
minant (5.3)  factorizes  into  similar  terms  with  different  k„  values.  Expanding 

(5.3)  in  terms  of  the  indices  m and  n shows  that  each  such  factor  is  given 


detj^^  (I-G^’V)  =[l-V(-l,-l)  G°(-l,-l)j  j^l-V(O.O)  G°(0,0)j  -i  ^G°  (0 , 0)G°  (-1 , -1 ; 


(5.4) 

In  this  expression  we  used  the  notation  det,  for  the  partial  determinant, 

is." 

with  a specific  k„. 

Substituting  the  explicit  expressions  for  the  surface  Green's  functions  G°(0,0) 
and  G°(-l,-l),  from  the  general  expression  (3.21),  into  Eq.  (5.4)  gives 


det.  (1-G  V) 

Kfi 


46  + M.  -H.(u)  -2U  ) 
la  a a 


-il(u)  -2U  ) 
a a a 


>FK- 


The  bound  state  energies  of  the  crystal  with  the  interface  are  given  by  the 

f (16-18) 

roots  of  the  equation  : 


det(l-G^V)  = 0 


(5.6) 


outside  the  shifted  bands  of  the  two  metals.  Using  the  explicit  expression 

(5.5)  we  see  that  this  is  exactly  the  same  condition  derived  earlier  from  the 
poles  of  the  interface  Green's  function,  Eq . (4.6). 


In  analogy  with  ordinary  scattering  theory,  one  defines  the  partial  phase 
shifts  by<'7-'8>: 


) 


■ 


n (E,jc„)  = Im  In  jdetk,,  (1-G°V) 


M-  + i (tiJ  -2U  ) 

a a 

+ i(oj  - 2 U ) 

a a a j 

(5.7) 


In  the  second  step,  we  applied  the  explicit  expression  (5.5)  and  used  the  iden- 
tity Iminf  = Argf,  where  Arg  denotes  the  argument  of  the  complex  function. 


If  we  take  the  determinant  of  (5.5)  with  respect  to  K„,  and  use  the  relation  (5.1) 
we  see  that  the  change  in  the  density  of  states,  due  to  the  formation  of  the 
interface,  is  given  by: 


A„ 
n llr 


Ap(E)  - § —^2  Ie  j Tl(E,k„ 


= £ 

n 


A„  , ^ l^a  - 2^3))  |^b^(-b-2U,)|)(5.8) 


In  this  expression  we  introduced  explicitly  the  g-fold  degeneracy  of  the  bands 
under  consideration.  We  also  note  that  A„  is  the  area  of  the  unit  cell  parallel 
to  the  interface  and  that  the  integration  of  k„  goes  over  the  two-dimensional 
Brillouin  zone  defined  by  the  crystals  structure  parallel  to  the  interface. 


Equation  (5.8)  can  be  used  to  derive  a relatively  simple  expression  for  the 
change  in  the  integrated  density  of  states.  This  is  given  by; 


AN(E)  = j Ap(E)  dE 


= £ 
T7 


A„ 


(2n) 


( 4B^+ 

LI  + i (ou  - 

a a 

■2^3)  k + 

-2U^) 

1 

i(uj  -2U  ) 
a a 

^b " 

i(c.b-2U^) 

Using  Eq,  (5.8),  and  integrating  by  parts,  it  is  not  hard  to  show  that  the  sin- 
gle-particle contribution  to  the  interface  energy,  (i.e.  the  change  in 

the  total  single-particle  energy  due  to  creation  of  the  interface  at  T=0), is 
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given  by  the  expression 


EAp(E)  dE 


•)  J J l|^  ) 


3 h 

where  E,  is  the  common  Fermi  energy  of  the  two  metals,  i.e.  E^  = ^ (E^  +E  ) and 

f t t g 

uj  and  a),  include  the  corresponding  shifts  in  the  energies,  according  to  Eq. 
a b 

(3,22). 

In  a similar  way,  one  can  determine  the  single-particle  contribution  to  the  change 
in  the  electronic  specific  heat.  This  change  is  given  by; 


= f EAd(E)  II  dE 


where  f(E)  is  the  Fermi-Dirac  distribution  function.  For  temperatures  much  lower 
than  the  Fermi  temperature  the  change  in  the  electronic  specific  heat  is  linear 
in  the  temperature  and  proportional  to  the  change  in  the  density  of  states  at 
the  Fermi  level  (assuming  no  Van  Hove  singularity  occuring  at  the  Fermi  energy). 
Explicitly  we  have 

_2 


where  K is  the  Boltzmann  constant.  For  the  interface  system  the  constant  y is 

D 

given  by 


7 - A„  |dk„  Arg 


a , ^ 4 a a M»,-2U,.|)Ci.ll) 


where  the  integrand  has  to  be  evaluated  at  the  common  Fermi  level. 
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APPENDIX  A 


In  this  appendix  we  derive  the  expressions  for  the  general  matrix  elements 

of  the  surface  Green's  function.  The  technique  we  apply  is  the  same  as  that 

(14) 

used  by  Kalkstein  and  Soven  to  derive  the  corresponding  diagonal  matrix 
elements . 

To  form  to  metal-vacuum  interface  we  start  from  the  bulk  metal  and  break  the 
bonds  between  the  two  adjacent  planes  n = 0 and  n = -1  . In  addition,  we 
change  the  self-consistent  potential  of  the  electrons  near  the  surface.  The 
perturbation  described  above,  which  is  necessary  to  create  the  surface,  has 
the  following  non-vanishing  matrix  elements  in  the  Uloch-Wannier  representa- 
tion : 


V(0,0)  = V (-1,-1)  = u 

(A.l) 

1 fl 

V(-1,0)  = V(0,-D*  = -Te 

(A. 2) 

With  this  perturbation  potential,  the  Dyson's  equation  for  the  surface  Green's 
function  takes  the  following  form: 

G(m,n)  = G°(m-n)  +G°(m+1)  V (-1 ,0)G  (0,n)-tG°  (m+l)V  (-1 , -1  )G  (-1  ,n) 

■♦G°(m)V(0,-l)G(-l,n)4G°(m)V(0,0)G(0,n)  (A.  3) 

In  writing  down  the  above  equation  we  used  the  property  that  due  to  transla- 
tional invariance,  the  bulk  Green's  function  G°  depends  on  the  indices  m and 
n only  through  the  different  m-n.  Taking  into  account  the  fact  that  the 
surface  Green's  function  G(m,n)  vanishes  if  m and  n refer  to  planes  on 
opposite  sides  of  the  interface,  we  can  solve  Eq . (A. 3)  for  the  general  matrix 
elements  of  the  surface  Green's  function.  This  is  given  by: 


A-1 


1 


(;(m,n)  = (.°(m-n)  + [g°  (m+l)V(-l  ,0)  + G°(m)V(0,0)|  G°(-n)  jl-G°(l)  x 

V(-1,0)-G°(0)V(0,0)|"^  (A. 4a) 


for  itijn-O 


G(m,n)=G°(m-n)+  |g° (m)V(0,-l)  + G° <m+l)V (-1 , -1 ) G°(-l  -n)  1-G°(-1)  x 


V(O.-l)  -cf(.l)V(-l,-l) 


-1 


(A4.b) 


for  m,ii  - -1 


G(m,n)  = 0 otherwise. 


(A, 4c) 


The  bulk  Green's  function,  G°(n),  was  derived  by  Kalkstein  and  Soven^^^^  and 


is  given  by; 


/ Oj-Hg,  { 

(— ) 


-in0 


Using  this  expression  and  the  expressions  (A.l)  - (A. 2)  for  the  perturbation, 
we  obtain  from  Eqs.  (A, 4)  that  the  matrix  elements  of  the  surface  Green's 
function  are  given  by: 


K».n) 


m+n  i|j,+(U)-2U) 


ig-(uJ-2U)  ) 


(A. 5a) 


for  mn,.'0 


!ml-f-!nl-2 


ig,+(q)-2U) 

iM.-(iu-2U) 


I (A. 


5b) 


for  m,n*.-l,  and 


G(m,n)  = 0 
otherwise . 


(A. 5c) 


A-2 


